Helioseismology 


In the 1960’s, it was found that turbulence in the Sun’s convective 
zone created waves below the surface. (Essentially, the convective 
motions were ringing a bell). Certain frequencies become amplified 
by constructive interference, leading to resonances, while others are 
destroyed by destructive interference. ‘The interpretation is that 
standing waves are trapped in cavities, with one end of the cavity 
being close to the stellar photosphere, and the other at a depth 
depending on the frequency. These waves manifest themselves as 
^ 1 m/s pulsations motions across the surface of the Sun. The 
typical period of these pulsations is ~ 5 min, though the power 
spectrum is complex, due to the beating of many, many periods. 


Consider two waves leaving A and arriving in B in phase. An os- 
cillation will be seen provided the paths differ by an exact number 
of wavelengths. In that case, the waves will reinforce each other; 
otherwise, they’ll interfere and have negligible amplitude. 


There are three types of oscillations that may occur in a star. 


p-modes 


p-mode oscillations are simply sound waves moving through gas, 
i.e., their restoring force is the local pressure. The speed of sound 
plus the ideal gas law implies 
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So the deeper into the star the wave goes, the larger the sound 
speed. In other words, for a wave traveling at an angle, the deeper 
wave front will always be traveling faster than the trailing wave 
front. This leads to refraction. 
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For a p-mode to create a standing wave, it must reflect off the 
stellar surface. If the pressure scale height, Ap, is larger than the 
wavelength of a mode, the wave has no trouble propagating and 
is transmitted. However, if the Ap is shorter than the wave, the 
gas readjusts on a timescale that is shorter than the wave’s period, 
and the wave is reflected. The upshot is that there is a maximum 


frequency to observable p-modes: 
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Since the star is in approximate hydrostatic equilibrium, and we are 
concerned with the stellar surface 
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p-modes are likely excited through stochastic acoustic noise from 
fluctuating turbulent pressure and gas pressure in the Sun’s convec- 
tive layer. Radiative losses, viscosity, and non-linear interactions 
between modes all work to dampen the oscillations. For low £ val- 
ues, a particular p-mode will decay exponentially on a timescale of 
days, suggesting 


1) whatever excites the p-modes happens on a timescale of a few 
days. 


2) destructive interference has plenty of time to destroy all the non- 
resonant frequencies, changing random convective noise into a very 
rich power spectrum. 


In general, low frequencies waves are damped out more slowly than 
high frequency waves. The lowest frequency modes may take months 
to decay. 


g-modes 


e-modes oscillations are driven by the negative buoyancy of verti- 
cally displayed material. Recall that the condition for convective 
stability is 
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This can be rewritten as 
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The buoyancy acceleration associated with the density perturbation 
is therefore 
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If we neglect the term associated with the chemical gradient, then an 
element which is displaced by a distance Ar in a region of stability, 
will move with Ar = Arge*"* with a frequency 
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This is the Brunt-Vaisala frequency, or buoyancy frequency. ‘The 
implication is that in regions of stability, material can oscillate and 
propagate g-modes, but in convective regions, displacements grow 
exponentially and no oscillations occur. 


In the Sun, g-modes can only occur in the central (radiative) regions. 
They are most likely excited by nuclear burning instabilities, and 
have periods greater than 50 min. 


f-modes 


f-modes are low-frequency surface waves propagating in their funda- 
mental (n = 0) mode. As with stars in the instability strip, simple 
gravity is the restoring force. They are similar to deep sea waves, as 
the medium acts as an incompressible fluid. Since these waves never 
penetrate more than ۸/4 into a star (where R is the stellar radius), 
they do not provide much information about the stellar interior. 


Pulsation Modes 


The pulsation modes of stars can be described by spherical harmon- 
ics and three integer numbers: 
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e The harmonic degree, £, is the number of nodes lines on the sur- 
face. Alternatively, one can think of / as the total number of planes 
slicing through the Sun. A purely radial pulsation has £ = 0. 


e the azimuthal order, m, gives how many of the surface node lines 
cross the equator. In effect, it describes how many planes slice 
through the Sun longitudinally. m ranges from —4 to +£. 





Examples of | = 3 models seen from different inclinations: 30° (top 
row), 60° (middle row) and 90° (bottom row). Blue represents areas 
of the surface moving inward; red shows regions moving outward, 
and white displays the nodal surface lines. From left to right, the 
columns display |m| = 3,2, 1, and 0. 
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Note: being acoustic waves, p-mode frequencies increase with an 
increase of the number of radial nodes (n) and degree £. But because 
buoyancy drives g-modes, their frequences decrease with larger n 
(but increase with £). 


The Sun 


The Sun has a convective envelope and a radiative core. As a result, 
the Sun’s g-modes are confined to the interior and are not directly 
observable. But ~ 10° separate p-modes and f-modes have been 
observed on the Sun’s surface. 
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lo measure p-modes in the Sun, one needs a long time series of 
observations (at least a month, ideally over several years) and obtain 
Doppler velocities (to better than ہ‎ 0.1 m/s) over the Sun's surface. 
(Because we are measuring radial oscillations, the data are much 
more easily obtained for the center of the Sun's disk that it's edge.) 


Typically, measurements are made once a minute. 


Velocity data are decomposed into spatial spherical harmonics, 


a = D —aegm(t)Y7" (6,9) (32.08) 
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where 0 is the co-latitude and © is the longitude. The aem am- 
plitudes contain all the information on the (¢,m) modes; a Fourier 
transform of a; m(t) then gives the frequencies of oscillations of the 
modes, a, (v). The final result is a power spectrum of three vari- 
ables: v, £, and m. 
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At a given £, each peak (ridge) represents a mode with different n 


(i.e., nodes in the radial direction): the higher the frequency, the 
higher the value of n, and the higher the depth of penetration. 
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The pattern shows a relation between the frequency and the mode 


4 1 
V f — Av (n ل 9 ل‎ Q -+ 3 ل‎ En 0 (32.09) 


where o is a function of the frequency determined by the properties 
of the surface layers, and رہ‎ is a small correction which depends 
on the conditions in the stellar core. Av is the inverse of the sound 
travel time across the stellar diameter: 


Av =2 0 3 i (32.10) 
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Note that from the virial theorum and the ideal gas law 
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Because the propagation of p-modes depends on the sound speed 
(which, through the equation of state, depends on p(r), u(r), and 
T(r)), and because each mode penetrates to different layer of the 
oun, the 2-D power spectrum can be used to obtain measure the 
internal structure and composition of the Sun. Moreover, 


e For a static Sun, the modes shouldn't depend on the direction 
of wave propatation. However, rotation breaks this degeneracy, as 


sounds will propagate faster in the direction of rotation than against 
the rotation. As a result, modes of the same (n,£) but different m 
will be split into 27 4- 1 frequency components. Specifically, a wave 
moving in the direction of rotation will have its frequency Doppler 


shifted by 
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where دن‎ is the rotation velocity and v, is the wave’s phase velocity. 
Similarly, a wave moving against the rotation will be shifted by 
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Hence the two waves will be split in frequency by 2Av. The result 
is that the entire standing wave pattern drifts slowly across the Sun 
in the direction opposite the solar rotation. By examining how this 
shift changes with (n,@), one can measure the Sun's rotation over 
90% of its mass. 





e For most applications, m splitting is measured, then removed, 
and the power spectra are averaged over the m states to produce 
the £ — v diagram. 


e The sound speed can be affected by localized phenomenon, such 
as the sunspots. (Their lower temperatures lower the sound speed.) 


Thus, one can “observe” activity on the far side of the Sun via the 
p- and f-modes, measure magnetic field strengths, etc. 


e Decause the sound speed depends on the mean molecular weight, 
helioseismology can be used to determine the Sun’s age: the larger 
the age, the greater the helium abundance, the larger the mean 
molecular weight, and the slower the sound speed. 


e Another way of looking at pulsations is through the use of two 
critical frequencies, the Lamb frequency and the Brunt- Vaisala fre- 
quency. ‘hey can be written as 
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Oscillations with w > Lg and w > N have pressure as the restoring 
force and are p-modes; on the other hand, if w < Ly and w > N, 
buoyancy is the restorer, and it is a g-mode. In stars like the Sun, 
the Lamb frequency decreases monotonically towards the surface, 
while N has a sharp peak near the center, and then rapidly drops 
to zero. 


e The solar oscillations, by definition, very slightly distort the shape 
of the Sun. 





Asteroseismology 


The pressure variations connected with the sound waves lead to 
fluctuations in temperature. This can change both the strength 
of the spectral lines being observed, and (very slightly) the total 
broadband intensity. Essentially, because of its p-mode oscillations, 
the Sun is a variable star which pulsates with an amplitude of ~ 
0.01 millimag. 
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Comparison of solar power spectrum from spectroscopy (red) and 
photometry (black). 


The changes in intensity open up the possibility of performing as- 
teroseismology on stars. ‘here are several types of stars for which 
such an analysis is viable; the most common of which are white 
dwarfs, pre-white dwarfs, solar-type stars, and subgiant stars. 
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White Dwarfs and Pre-White Dwarfs 


Aside from the Sun, the first non-radial pulsators to be analyzed 
by asteroseismological techniques were white dwarfs. Since white 
dwarfs (and planetary nebula nuclei) are hot, they do not have 
convective envelopes. This means that a) there is no driver for p- 
modes, and b) g-modes can propagate to the surface. As a result, 
white dwarfs (and pre-white dwarfs) that are in an instability strip 
(DAs for the partial ionization of helium, DBs for partial ionization 
of doubly-ionized helium and DOs and planetary nebulae nuclei 
for the partial ionization of carbon and oxygen) can be non-radial 
pulsators. 


DA white dwarf variables (DAVs) are called ZZ Ceti stars; DOVs 
are generally called GW Vir or PG 1159 stars. DB variables are 
usually just called DBVs. 
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White dwarfs can pulsate with amplitudes of ~ 0.05 mag, and much 
of their power coming out in periods of ~ 10 min. (For pre-white 
dwarf planetary nebula central stars, this period is about ~ 25 min.) 
Moreover, often times, these objects have two principle modes that 
are moderately close together. That means that the periods beat: 
every few hours the positive and negative excursions add up, cre- 
ating variations of almost ~ 0.2 mag. So detecting multi-mode 
pulsator from the ground is easy! However, in order to extract in- 
formation about the stars, one needs to resolve all the modes. 


A rule of thumb: to resolve two frequencies separated by Av, one 
needs to observe (uninterrupted) for 17/7 seconds. So to identify two 
modes separated by 2 uL Hz, one needs a month or so of continuous 
observations! 


In 1990, the Whole Earth Telescope was formed from a network of 
small telescopes around the globe. Since then, other networks of 
telescopes have investigated a small handful of objects. 
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The light curve of white dwarf PG 1159-035 over the central ~ 6 
days of a month-long campaign with the Whole Earth Telescope. 
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The power spectrum of PG 1159-035. 
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Recall that g-modes propagate at the Brunt-Vaisala, or buoyancy, 
frequency which, in the absense of a chemical gradient (see 32.07) 
is 


Ne = geo [Vad — Viad] (32.16) 


Unlike p-modes, which exhibit à periodicity in frequency, g-mode 
pulsations are more-or-less periodic in their periods. In the limit of 
high-n, the asymptotic form of the period relation is 
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So the period spacing for £ = 1 modes will be a factor of 4/3 larger 
than for £ = 2. Ilo is related to the buoyancy frequency, N, by 
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where rı and rə are the boundary regions for the mode. For white 
dwarfs, this quantity depends primary on the total stellar mass, 


M ۱ ۱ 
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The mass of a well-observed white dwarf can therefore be measured 
to 3 significant digits. Other information obtainable include the 
rotation rate, the thickness of the non-degenerate surface layers, 
and the composition of both the main body of the star and the 
surface layer. 





Finally, there is also the possibility of observing actual stellar evo- 
lution via temporal changes in the periods of white dwarfs. As 
described above, g-modes have a period periodicity. If k is the wave 
number, then the buoyancy is related to period by 
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For roughly isothermal cores (as might be expected from a white 
dwarf with Urca process neutrino cooling), one can write down the 
approximate dependence of period on core temperature and stellar 


radius 
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dt 2 dt dt 
At first glance, this relation looks problematic, since, as a hot white 
dwarf evolves, it should both cool and contract, and those two ef- 
fects work in opposite directions. Fortunately, the contraction of 








(32.20) 


white dwarfs (and pre-white dwarfs) should be minor, implying that 
the g-mode periods should decrease with time. This makes empir- 
ical sense, as pulsating planetary nebulae have principle periods of 
~ 25 min, while white dwarf periods are 5 to 10 min. However, 
the attempts to directly observe white dwarf cooling have thus far 
produced ambiguous results. 
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Solar type Stars 


Over the years, there have been many attempts to detect the p- 
mode oscillations of nearby solar-type stars through ground-based 
observations. But (with one or two exceptions), it was not until Ke- 
pler that such observations became feasible. Now, with the Kepler 
and CoRoT satellites, we can analyze the p-modes of large num- 
bers of stars and derive their stellar parameters from their general 
oscillation properties. This is called “Ensemble Asteroseismology.” 


The power spectrum of most solar-type stars looks similar to that of 
the Sun. There is a regularity in the frequency separation of power 
related to the sound travel time across the stellar diameter 


Ay fü z) x (p)? (32.10) 


This is generally called the “large separation”. Since c, « رو‎ zi 


this periodicity generally scales with the square root of the mean 
density of star: as the star gets larger, Ar decreases. 


There is also another periodicity known as the “small separation" 
in the power spectrum peaks 











Óv, = Vn e — Vn—1,042 = (4£ aS 6) Do (32.21) 
where 
Av |c.(R) R de, dr 
Do = — — — 32.22 
à م4‎ ۸,۰۸۱ R | dr ( 





This small frequency separation is most sensitive to the sound-speed 
gradient at small r (ie., the core), which in turn is sensitive to 
the core's chemical composition gradient and hence to the star's 
evolutionary state. 


Finally, the overall pattern of power in the oscillations looks ap- 
proximately Gaussian, with the bulk of the power coming out at a 
frequency, Vmax. This frequency is likely related to the maximum 
frequency cutoff (32.04). If so, then as a star evolves towards the 
subgiant branch and its surface gravity decreases, Vmax Should move 
towards lower frequencies. 
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Power spectra for a wide binary system, 16 Cyg A (left) and 16 Cyg 
B (right). Note the Gaussian-like peak at high-frequency, and the 
periods in frequency space. 


























If the peak in the distribution of frequency modes is related to the 
cutoff frequency (32.04), and through (32.12) Av « (p)!/? then a 
simply scaling to the Sun yields 
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Putting these two equations together yields 
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For the Sun, Vmax ~ 3090 wHz, Av ~ 135.1 wHz, and Teg ~ 5777 K. 
So, without any complex modeling, measurements of Vmax and Av 
can produce estimates of M and R than are better than ~ 10%. 


We can also use the large and small separations in the power spec- 
trum peaks to estimate the age of a solar-type star. Such a mea- 
surement requires a set of evolutionary models (to relate quantities 
such as the mean density and core mean molecular weight to age), 
but are fairly straightforward. Below is one such set of models for 


solar-metallicity stars. 
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Diagram showing large separation versus small separation for a 
series of evolutionary models calculated for different masses, but 
fixed metallicity. Solid lines are evolutionary tracks for decreasing 
core’s hydrogen abundance. Dashed black lines indicate models with 
same core’s hydrogen content decreasing from X, = 0.7 (ZAMS) to 
X. = 0.1. The Sun and other well-observed stars are marked. 






































Note: one complication in interpreting the power spectrum of sub- 
giants and red giants is “mixed mode” oscillations. In solar-type 
objects, interior g-modes have low frequency, while most p-modes 
have a relatively high frequency, so the two modes do not interact. 
However, as stars evolve, the large chemical gradients that are built 
up and the large density gradient between the core and the envelope 
can increase the buoyancy frequency. ‘These interior g-modes can 
then interact with p-modes of roughly the same frequency, creating 
modes of mixed character. As a result, frequencies become shifted 
from the regular Av spacing, and the star’s power spectrum has a 
quite complicated appearance (at least in the region of g-mode and 
p-mode frequency overlap). 


Most p-mode-based estimates of stellar mass, radius, and age are 
based on ensemble asteroseismology. In fact, ensemble asteroseis- 
mology in open clusters are particularly useful, as the technique 
naturally works best when the stars being considered are similar. 
For example, by determining the masses of red giants as a func- 
tion of luminosity in a single cluster, one can measure mass-loss 
directly, and not be affected by systematic uncertainties due to stel- 
lar metallicity, and the like. (In theory, one could have also done 
this with globular clusters, if there were a suitable such object in a 
Kepler field. But crowding and distance more-or-less excluded this 
possibility). 


An alternative and more precise way of determining stellar param- 
eters is through the analysis of individual oscillation frequecies, 
rather than the globular parameters of Vmax, Av, and di. ‘This 
requires the computation of detailed stellar models for each star, 
in order to compare theoretical oscillation frequencies with the ob- 
served frequencies. Obviously, this involves a lot more work, but it 
typically doubles the precision of a mass, radius, and age estimate. 


The comparison of observed frequencies versus theoretically com- 
puted frequencies is often performed using an “échelle” diagram. 
For the diagram, a power spectrum is divided into segments of 
length Av and each segment is plotted right above the other. If 
the Av spacing were perfectly regular, the result would be a ver- 
ticle line in the diagram. Curvature or other departures from this 
line contain the detailed information about the model. (In particu- 
lar, points which are well off the ridge lines represent mixed-mode 
frequencies. Since mode-mixing increases as the star evolves, the 
presence of such oscillations provides information about the evolu- 
tionary state of the star. In fact, by matching the observed and 
theoretical frequencies in an échelle diagram, it is possible to fix the 
mass and radius of a star to ~ 2% and the age to ~ 7%. 
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An échelle diagram for the Sun (left) and red giant star KIC4351319 
(right). The filled points are observations, while the open shapes are 
theoretical calculations. Points will large deviations from the ridge 
lines represente mixed-mode oscillations. Circles represent modes 
with £ = 0, triangles show £ = 1, squares £ = 2, and diamonds £ = 3. 
The distance between two adjacent frequency columns represents 
the small separation, ۰ 








Precision of Parameter Estimates 


Radius Mass Age 
Sun-Like Main Sequence 


CMD Diagram - 0 ~ 15% 
Ensemble Asteroseismology ~ 6% ~ 6% دا‎ 
Individual Frequencies 0.3% 0.8% 3.7% 


Evolved Stars 


CMD Diagram oe - 
Ensemble Asteroseismology ~ 1596 ~ 15% see 
Individual Frequencies 1% 2% 7% 


Stellar Rotation 


For a static star, the orientation of the star is irrelevant, and all 
the m modes are degenerate. But, just like the Zeeman and Stark 
splitting in atoms, the inclusion of a preferred direction (in this case, 
the rotation axis) causes the m modes to be split into 2/4- 1 different 
components, 1.e., 

Vn, £m = Vn lym — Vn,£0 (32.27) 


As long as the rotation is slow enough so that the centrifugal force is 
negligible, (32.13) and (32.14) hold. With a long enough baseline (as 
provided by Kepler), the splitting of these modes can be measured, 
and, by comparison to predictions of stellar models, a star’s rotation, 
as a function of radius, can be determined. This is best done with 
sub-giants and red giants because 


1) Mode-mixing allows us to constrain the star’s interior g-modes. 


2) The evolution off the main sequence (and the development of 
shell burning) decouples the core from the surface, allowing both 
regions to evolve independently. Due to core contraction and en- 
velope expansion, the core rotation should speed up, the surface 
rotation should slow down, with the hydrogen-burning shell divid- 
ing the two zones. ‘This is confirmed by observations, though the 
theoretical models predict rotation rates that are an order of mag- 
nitude faster than what is inferred from seismology. 
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lwo models for the internal rotation of RGB star KIC 4448777 
based on asteroseismology. ‘he dashed line shows the hydrogen- 
burning shell, while the shaded region indicators regions where no 
solutions could be found. In general, the core of RGB stars rotates 
8 to 17 faster than the outer layers, and the mean rotation slows as 
the star ascends the giant branch. 








Pulsation Amplitudes 


There is still no theory to predict the amplitude of stellar pulsa- 
tions. (While there are arguments for how the excitation of pulsa- 
tions scales with stellar conditions, our understanding of damping is 
poor.) Using empirical Kepler data, one can express the oscillation 
amplitude as a function of luminosity and mass, i.e., 
Le 
7ك‎ 
with s ~ 0.9 and t ~ 1.9. Alternatively, the pulsation amplitude can 
be expressed using a “flicker” parameter which scales with surface 


gravity. 
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Either way, it appears that as stars evolve off the main sequence 
towards the giant branch, their oscillation amplitudes increase. 





